NATICK7  DEVELOPMENT  CENTER. 
NATICK,  MASSACHUSETTS  01760 


iCHNICAL  REPORT 
75‘99-AiflEL 


THEORETICAL  AND  EXPERIMENTAL 
INVESTIGATION  OF  THE  BUCKLING  OF 
RECTANGULAR  CYLINDERS 


DDC 

Ji  rEB  24  1976 


February  1976 


Approved  for  public  releatej 
distribution  unlimited. 


Aero-Mechaucal  Engineering  Laboratory 


iimmk, 


s a iiv 

1 s I 

M ■ 

Sra  o 1 

I 1 

^ ? I 

■ .ji ' , jhkw’o 

Bg  < 1 

I r\ 

SI 

\*\  si 

W$£ 

f OM 

for  ^Hc 

Citation  of  «“ 

constitute  an  official  mu 

iMO  of  such  items* 

,*  ohm  no  longor  n«a*».  “»* 

Ttantrov  this  report  ufaen  no  ion* 

^tunTit  to  the  originator. 


i 


li 


:■ 


#ASfefea^M4Ml« 


^rm4-' 


yapfE 


IM 

RC 

l-  '’■' 
K&' 


iffi 


r^v.1 


iti 

®SSS 


gEfSg 

5*.  >.0*  SJfcjt 


fjKVj 


1 1 


sp; 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whon  Dmto  Entrrrdj 


RE  PORI  DOCUMENTATION  PAGE 

r REPORT  NUMBER  (/V "7  '*  r'r"J' 

Techoicaf-fteport)75-99-AMEL  / 


2.  govt  AvcEssroimoq 

LlI 


l 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


j Rwmwrsrmrss  NUMBER 

f\ 


t crktnTp 


TJTI  E (onri  SubUtUt- 


JHEORETICAL  AND  EXPERIMENTAL  INVESTIGATION 
OF  THE  BUCKLING  OF  RECTANGULAR  CYLINDERS, 


7.  AUTHORfF) 


\j  Earl  C.|Steeves  J 


s.  PERFORMING  ORG.  REPORT  NUMBER 
8.  CONTRACT  OR  GRANT  NUMBERf.) 


t.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

U.S.  Army  Natick  Developmpnt  Center  1 
Aero-Mechanical  Engineering  Laboratory 
Natick,  Massachusetts  017U0 


ln.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  8 WORK  UNIT  NUMBERS 


6.2;1J6627080552;05;017 


II.  CONTROLLING  OFFIC-  NAME  AND  ADDRESS 

U.S.  Army  Natick  Development  Certer 
Aero-Mechanical  Engineering  Laboratory 
Natick,  Massachusetts  01760 


ujStpjQfli.cua*  - 

FebnaMW76 


j'11  UMMBBn  >r  rnnn 

38 


l 


U MONITORING  AGENCY  NAME  A ADDRESS^!/  ditto  rwi<  from  Controlling  Olllct) 


15.  SECURITY  CLASS,  (ol  thto  roport) 


Unclassified 


15  DISTRIBUTION  STATEMENT  (oi  thldMopoTIJ, 


V / 


~ ,j_ 


is.,  declassification/downgrading 
schedule 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  <o  / tho  oiolroctonlorrTlntVoclr  16,  If  from.  Raport) 

Sr  - 1-7-  L ■ ■ ■ : -e'-  <T. 5 / 

Sr-  i-  L - ■ '■  * r - 7 • / 


i yt 


' / 


- a 


18.  SUPPLEMENTARY  NOTES  

/-J-  LL  - - 

/-  r - 


/ 


(7  t 


19.  KEY  WORDS  (Contlnum  on  row  tom  mi  dm  It  nacmaaary  and  Idmntlty  by  block  numbar) 


CYLINDERS 

RECTANGULAR  BODIES 
EQUATIONS 
BUCKLING 
GEOMETRY 


LOAD  CONTROL 
LOADING 

STRESS  ANALYSIS 

CONTAINERS 

PACKAGING 


CONSTANTS 

RATIOS 

EIGENVALUES 

MATRICES 

VECTORS  (MATHEMATICS) 


20.  ABSTRACT  (Contlnum  on  rmvmrma  ml  dm  tt  nmcmmmary  and  Idmntlty  by  bio  .k  numbar) 

^ Classical  plate  buckling  theory  is  used  to  n odel  the  buckling  problem  for  rectangular 
cylinders.  Solutions  to  the  resulting  equations  are  found  directly  for  the  case  of  a uniformly 
loaded  cylinder  and  by  approximate  techniques  for  the  nonuniformly  loaded  cylinder.  Results 
from  a series  of  experiments  to  determine  the  buckling  load  are  presented  and  compared  with 
predictions  from  the  theory.  These  comparisons  reveal  the  need  to  examine  geometrically 
nonlinear  buckling  problems  to  mo>  j closelv  nodel  the  situatir  i present  during  the  compression 
loading  of  rectangular  cylinders.  In  additit  • the  need  to  include  inelastic  behavior  for  cylinders 
^ofJow<>as£ectj|atioiJSijJso>jjeveale^ 


DD  1473  EDITION  OF  I NOV  *»  IS  OBS^SeTE 


Unclassified 


SECURITY  CLASJ'FrCATION  OF  THIS  PAGE  fWh»o  Dot « Bntorod) 


/ 


y 


i 


PREFACE 


The  results  presented  in  this  paper  were  developed  during  an  investigation  of  the 
behavior  of  package  and  container  structures  under  load.  The  wcrk  was  part  of  a program 
to  formulate  a computer-aided  design  procedure  for  packaging  systems  under  Task 
/1J6627080552  05/017,  Establishment  of  Design  Criteria  for  Containers,  and 
^1E662703A090-04  Design,  Analysis  and  Optimization  of  Struc  ures.  This  last  task  was 
sponsored  by  the  AMC  Computer-Aided  Design  and  Engineering  Program. 
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THEORETICAL  ANk.  'PERIMENTAL  INVESTIGATION  OF 
THE  BUCKLING  OF  RECTANGULAR  CYLINDERS 


INTRODUCTION 

A review  of  the  package  design  problem  reveals  a problem  involving  many  variables 
and  a great  quantity  of  data  restive  to  these  variables  and  the  design  parameters.  The 
complete  design  problem  involves  the  interaction  of  such  parameters  as  shipping  route, 
its  associated  hazard  environment,  and  characteristics  of  the  item  being  packaged,  such 
as  fragility  and  cost,  with  the  package  and  container  design  parameters.  A problem  of 
such  complexity  is  well  suited  to  the  computer-aided  design  approach.  One  of  the 
requirements  for  the  application  of  this  approach  is  the  ability  to  specify  a package  or 
container  design  which  will  protect  the  contents  against  the  expected  shipping,  handling, 
and  storage  environment.  One  element  of  this  environment  is  the  load  imposed  on  the 
package  or  container  and  the  work  reported  here  is  an  investigation  of  the  behavior  of 
recta. igular  packages  and  containers  under  one  of  these  loading  conditions,  namely,  the 
compressive  loading  due  to  stacking  during  shipment  or  storage. 

The  stacking  problem  is  treated  as  the  buckling  problem  for  rectangular  cylinders 
subject  to  axial  compression,  the  objective  being  to  develop  a program  which  would  allow 
specification  of  a design,  given  the  expected  stacking  loads.  To  accomplish  this  objective, 
a procedure  to  compute  the  buckling  loads  of  rectangular  cylinders  was  developed  using 
the  theory  of  structural  mechanics.  In  addition,  tests  were  conducted  to  experimentally 
determine  the  buckling  load  of  a series  of  rectangular  fiberboard  cylinders.  Comparison 
of  these  experimental  results  with  the  computed  values  provided  a measure  of  the  usefulness 
of  the  theory.  Although  the  experimental  work  reported  here  deals  solely  with  fiberboard 
cylinders,  the  theory  is  applicable  to  any  orthotropic  material. 

Previous  work  related  to  this  problem  has  largely  been  directed  toward  developing 
empirical  correlation  formulas  between  the  container  compressive  strength  and  such 
properties  as  ring  crush  strength,  basis  weight  and  Young's  Modulus.  The  results  presented 
in  reference  1 typify  these  efforts  in  which  some  of  the  notions  of  structural  mechanics 
are  utilized  but  a direct  attack  on  the  problem  is  not  made.  The  nost  direct  attack 
on  the  problem  known  to  the  author  is  reported  in  reference  2.  In  this  work  the  theory 
for  the  buckling  of  flat  plates  war  utilized  to  predict  the  strength  of  rectangular  containers; 
however,  no  attempt  was  made  to  treat  an  entire  package. 


'P.E.  Shick,  and  N.C.S.  Chan;  Top-to-Bottom  Compression  for  Double  Wall  Corrugated 
Boxes;  TAPP  I,  Vol  48  No  7;  1965;  pp  423-430. 

2B.S  Angell,  and  P.R.  Paslay;  Determination  of  the  Stacking  Strength  of  Corrugated 
Fiberboard  Containers;  Proceedings  of  Society  of  Experimental  Stress  Analysis,  Vol  XVI, 
No  1 pp  109-116. 


A.t-7  £ 


“-‘■OC.HIQJ, 


pPTt  . 


5 


THEORETICAL  PREDICTION  OF  THE  BUCKLING  LOAD 


The  problem  to  be  considered  here  is  that  of  the  strength  of  a rectangular  container 
under  stacking  loads.  The  main  load-carrying  elements  of  the  container  are  the  four  sides 
which  support  the  load  in  compression.  Because  cf  the  compressive  load  the  buckling 
mode  of  failure  must  be  considered.  To  simplify  the  theoretical  treatment  of  this  mode 
of  failure  the  stacking  problem  is  represented  by  c rectangular  cylinder  subjected  to  axial 
compression  as  illustrated  in  Figure  1.  In  this  section  the  theoretical  aspects  and  solutions 
or  tvo  approaches  to  the  problem  will  be  presented.  The  first  approach  is  a direct 
solution  of  the  governing  equations  leading  to  a transcendental  eigenvalue  problem  for 
the  buckling  load.  The  second  approach  is  a finite  difference  solution  which  leads  to 
an  algebraic  eigenvalue  problem. 

Direct  Solution  of  Governing  Equations 

Governing  Equation  and  Solution: 

The  formulation  of  the  problem  of  the  buckling  of  rectangular  cylinders  presents 
some  mathematical  difficulties  because  of  the  discontinuity  in  the  tanget  to  the  cylinder 
as  one  proceeds  around  its  perimeter.  This  difficulty  is  overcome  by  treating  each  of 
the  four  sides  of  the  cylinder  as  flat  plates,  as  shown  in  Figure  2,  and  accounting  for 
the  interaction  of  the  sides  through  boundary  conditions  expressing  continuity  of 
displacement  and  equilibrium  of  moments  at  the  corners.  To  accomplish  this  a coordinate 
system  is  established  on  each  of  the  sides  as  shown  in  Figure  2.  The  equations  governing 
the  buckling  of  each  of  the  plates  subject  to  an  axial  compressive  stress  Nx  are: 


33 4w;  34w:  34w:  32w: 

D, , + 2(D,  j + 2D33)  + D22  Nx  - 0 

3x4  3x23y?  3y4  3x2 

(1) 

A complete  development  of  these  equations,  which  are  written  for  an  orthotropic  plate 
whose  principal  directions  are  parallel  with  the  Xi  ,Yi  directions,  is  given  in  reference  3. 
In  equation  (1)  i takes  on  the  values  1-4  giving  an  equation  for  each  plate.  This  equation 
is  nondimensionalized  by  dividing  by  DM  and  letting  x=a£  and  yj=arjj: 

34Wj  34w;  34  W;  32w; 

+ 2e, + e, + N - 0 

3£4  3£23t?|  3t/!  3$2 

(2a) 

D j 2 .D  3 3 

e,  « — + 2— 

Dm  Dii 

(2b) 

D2  j 

- = D„ 

(2c) 

N - N*"’ 
D,, 

(2d) 

3R.F.S.  Hearmon;  An  Introduction  To  Applied  Anisotropic  Elasticity;  Oxford  University 

Press;  1961. 
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oads  on  rectangular  container 


Figure  2 Coordinate  Systems  used  in  Analysis  of  Rectangular  Cylinder 


In  writing  equation  (2)  account  has  been  taken  of  the  compressive  nature  of  Nx  by  changing 
the  sign  of  the  term  containing  this  parameter.  It  is  assumed  that  ‘he  boundaries  £~0 
and  {=1  are  simply  supported,  that  is,  the  displacement  and  the  momant  vanish  at  these 
boundaries.  With  this  assumption  the  normal  displacements  w;  ^re  expressed  as: 


Wj(£,r?j)  = Fj(rjj)  Sin  mrf 

and  equation  (2)  is  reduced  to 

, dJ  Fj  d4F; 

(nir)4  Fj  -2e!  (n;r)J + e2 - N(njr)JFj  = 0 

drj?  dr?f 

The  solution  of  this  equation  is 


(3) 


(4) 


Fj(tjj)  = AjSinhorTjj  + BjCoshorrjj  + CjSinflrjj  + DjCos/fyj  (5) 

where 

a = 

0 = 

This  gives  the 

w;  = Sin(nrr<)  [AjSinh(arjj)  + BjCoshfcn?-,)  + CjSini^rjj)  + DjCos(0rjj)J  (?) 

With  the  constants  of  integration  and  the  buckling  eigenvalue  to  be  determined  by 
application  of  the  conditions  joining  the  corners. 

Joining  Constraints: 

The  developments  up  to  this  point  have  considered  the  sides  of  the  rectangular  cylinder 
as  though  they  were  not  joined.  In  order  to  complete  the  solution  it  is  necessary  to 
state  the  conditions  joining  the  sides  and  apply  these  conditions  to  the  solution  obtained 
above. 


(6a) 


(6b) 


:omplete  solution  as 


Because  of  the  relatively  greater  inplane  stiffness  of  the  plates  it  is  assumed  that 
at  each  corner  the  normal  displacement  vanishes;  this  is  expressed  mathematically  as: 
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(3) 


w,(£,  b/a)  = 0 
w,  (£,  -b/a)  = 0 
w2(£,  c/a)  = 0 
w2(£,  -c/a)  = 0 


w3(£,  b/a)  = P 
w3(£,  -b/a)  - 0 
w4{£,  c/a)  = 0 
w4(£,  -c/a)  = 0 


Equations  (8)  must  be  satisfied  for  all  £. 

application  of  these  conditions  to  the  solution  allows  the  determination  of  eight  of  the 
sixteen  constants  in  equation  (7)  yielding  the  following  for  the  solution: 

w,  = C,  [Sin(0i?, ) + f2Sinh(aij, )]  + D,[Cos(0t?,)  + f,  Cosh(aj?, )] 

Sin(njr£) 


w2  = C2[Sin(0tj2)  + f4Sinh(orT72 ) ] + D2[Cos(0t?2)  + f3Cosh(ai?2)] 
Sin(ntr£) 


(9) 


w3  = C3  [Sin(/3i?3 > + f2Sinh(ccr?3)]  + D3  [Cos(/3i73 ) + f,  Cosh(ori73)] 

Sin(n7r£) 

w4  = C4  [Sin(/3i?4 ) + f4Sinh(ar?4)]  + D4  [Cos(0r}4)  + f3Cosh(ar?4)] 

Sin(n7r£) 

In  equation  (9)  the  constants  f,  have  the  following  definition: 
f,  = -Cos(0b/a)/Cosh|ab/d) 
f2  = -Sin(0b/a)/Sinh(ab/a) 
fj  = -Cos(/3c/a)/Cosh(crc/a) 
f4  ■=  -Sin(0c/a)/Sinh(ac/a) 

The  remaining  eight  joining  conditions  to  be  applied  express  the  continuity  of  first  and 
second  derivatives  of  the  displacements  at  the  corners  as: 


(10) 


w,'(t  -b/a)  - w2'(£,  c/a)  * 0 
- w’,'(£,  -b/a)  + w2  (£,  c/a)  = 0 


(11) 
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W 2(£, 

-c/a) 

- w3'(£,  b/a) 

- w'2($, 

-c/a) 

+ w'3'(£,  b/a) 

W3(£, 

-b/a) 

- w'4({,  c/a) 

- w'a'U, 

-b/a) 

+ w4'(£,  c/a) 

wi(£, 

-c/a) 

- w,'i£,  b/a) 

- wi(£, 

-c/a) 

+ wi'(£  b/a) 

for  all  £. 


(ID 


Where  the  primes  denote  differentiation  with  respect  to  the  appropriate  variable. 
Substitution  of  equations  (9)  into  (11)  yields  the  following  set  of  homogenous  equations 
written  in  matrix  form: 

[U]  (cj  =0  (1?) 

In  equation  (12)  { c]  is  a vector  whose  elements  are  the  const  C, , D, , C2 , D2 , C,, 

Dj,  C4,  D4  and  the  nonvanishing  elements  of  the  matrix  [U] , v,:  : .h  is  of  order  are 

given  in  Appendix  I. 

A non  trivial  (cj  which  satisfies  equation  (12)  exists  only  if  the  matrix  U is  singular; 
thus  one  seeks  the  values  of  N which  cause  U to  be  singular  and  the  corresponding  solution. 
These  are  respectively  the  buckling  eigenvalue  and  buckling  mode  shape. 

Extraction  of  Buckling  Eigenvalues: 

As  indicated  =ibove,  the  values  of  N which  cause  the  matrix  U to  be  singular  are 
sought;  however,  the  parameter  N appears  tmscendentally  in  the  elements  of  U so  a 
direct  attack  is  not  possible  and  a search  and  interpolation  scheme  must  be  used.  The 
scheme  used  to  obtain  the  results  presented  in  this  report  is  basad  on  the  theorem  presented 
and  proved  in  reference  4 which  deals  with  the  solution  of  sets  of  homogenous  equations 
in  which  one  or  more  of  the  equations  are  redundant.  The  results  presented  in  reference  4 
are  motivated  by  the  situation  in  which  the  number  of  equations  is  less  than  the  number 
of  unknowns;  however,  the  extension  to  the  present  case  is  straightforward. 


4W.C.  Walton,  Jr.,  and  E.C.  Steeves;  A New  Matrix  Theorem  and  Its  Application  for 
Establishing  Independent  Coordinates  for  Complex  Dynamical  Systems  with  Constraints; 
NASA  TR-R-326;  1969. 
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According  to  the  theorem  proven  in  reference  4,  the  solution  of  a system  of 
homogenous  equations  like  (12)  can  be  expressed  as  a linear  combination  of  the 
eigenvectors  corresponding  to  the  eigenvalues  of  the  product  matrix  lP"U  having  the  value 
zero.  When  the  matrix  U is  square  and  nonsingular  U^U  will  have  no  zero  eigenvalues 
and  equation  (12)  will  have  no  non-trivial  solution.  Thus  determination  of  the  buckling 
eigenvalues  is  reduced  to  finding  the  veiue  of  N which  causes  one  of  the  eigenvalues  of 
l)Tu  to  vanish.  The  fact  that  the  singularity  of  U will  cause  one  of  the  eigenvalues 
of  llTu  to  vanish  results  from  the  fact  the  singularity  implies  a redundant  equation,  and 
it  is  shown  in  reference  4 that  a redundant  equation  manifests  itself  in  a vanishing 
•igen' alue  of  U^U.  It  should  be  noted  that  although  U is  square,  it  is  not  symmetric; 
because  of  this  a procedure  based  on  the  eigenvalues  of  U was  not  adopted  because  of 
the  question  of  realness  of  eigenvalues  that  arise  with  nonsymmetric  matrices  The  search 
aspect  of  the  solution  scheme  is  carried  out  by  computing  the  product  of  the  eigenvalues 
of  UTu  for  successive  values  of  N until  a value  of  N close  to  the  solution  is  found. 
Because  the  product  matrix  U^U  is  positive  semi-definite  the  eigenvalues  are  always  greater 
than  or  equal  to  zero.  Thus  the  product  of  the  eigenvalues  is  greater  than  or  equal  to 
zero  giving  the  behavior  of  the  product  of  the  eigenvalues  as  a function  of  N near  a 
solution  as  shown  in  Figure  2a.  The  usual  technique  of  looking  for  a change  in  sign 
to  detect  nearness  to  a solution  is  not  applicable  in  this  case  and  the  alternate  procedure 
of  looking  for  changes  in  the  sign  of  the  difference  of  successive  values  of  the  eigenvalue 
product  is  used.  This  is,  in  effect,  looking  for  changes  in  the  slope  of  the  function  shown 
in  Figure  2a.  Once  the  location  of  a root  is  found  through  this  search  procedure  an 
interpolation  procedure  is  used  to  find  the  value  of  the  root  to  within  a specified  degree 
of  accuracy.  A quadratic  interpolation  is  used  because  of  the  behavior  near  the  root. 
When  a root  has  been  located  three  pairs  of  values  (N,  P)  where  P is  the  product  of 
the  eigenvalues  of  U^U,  are  known  and  are  used  to  determine  a quadratic  function.  The 
next  trial  value  of  N is  that  value  which  causes  the  slope  of  this  quadratic  function  to 
vanish.  Using  this  value  of  N the  corresponding  value  of  P is  computed,  and  this  (N, 
P)  pair  is  used  in  conjunction  with  the  two  most  recent  pairs  to  obtain  a new  quadratic 
interpolation  function  which  yields  a new  trial-value  of  N.  This  interpolation  procedure 
is  continued  until  the  normalized  difference  Letween  two  successively  computed  values 
of  N is  less  thar  a specified  convergence  criterion. 


Two  questions  arise  regarding  the  use  of  this  technique  for  finding  the  roots.  The 
first  concerns  the  hehavior  when  the  matrix  U is  poorly  conditioned,  that  is,  when  the 
determinant  of  U approaches  the  value  zero  for  values  of  N not  near  a solution.  In 
such  a case 

Det.  U = f 

where 

e < 1 
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The  present  technique  utilizes  the  matrix  llTu  for  which  the  determinant  in  the  case 
of  a poorly  conditioned  matrix  is 

Det  UTU  2 e1 

Thus,  in  comparison  to  the  commonly  used  techniques  based  un  the  determinate  of  U, 
the  present  technique  is  more  prone  to  conditioning  problems  because  of  the  above 
quadratic  dependence  on  e.  This  problem  was  not  encountered  in  any  of  the  computations 
carried  out  in  this  work  indicating  that  the  equations  expressed  are  clearly  linearly 
independent.  The  second  question  concerns  the  possible  detection  of  false  solutions.  Such 
false  solutions  correspond  to  values  of  N for  which  the  slope  vanishes  with  a non  zero 
value  for  the  product  of  eigenvalues.  No  automatic  procedure  was  developed  to  rule 
out  this  possibility  but  in  each  case  the  eigenvalues  corresponding  to  a solution  were 
examined  to  verify  that  at  least  one  of  them  vanished.  Because  of  round-off  it  is  necessary 
to  examine  the  eigenvalues  as  opposed  to  their  product. 

Results  using  this  technique  are  presented  in  a subsequent  section  of  this  report. 
Finite  Difference  Solution 

The  direct  solution  of  governing  equations  described  in  the  preceding  section  is 
computationally  inefficient  because  of  the  indirect  approach  to  eigenvalue  extraction  and 
is  limited  to  the  treatment  of  uniformly  loaded  cylinders.  In  an  effort  to  improve 
computational  efficiency  and  develop  a method  capable  of  treating  nonuniformly  loaded 
rectangular  cylinders,  an  approximate  solution  based  on  finite  difference  procedures  was 
developed. 

Energy  Minimization  Procedure: 

Although  finite  difference  procedures  are  generally  applied  directly  to  differential 
equations,  it  is  possible  to  use  this  approximation  technique  to  reduce  an  integral  energy 
minimization  principle  to  a system  of  algebraic  equations  without  first  developing  the 
differential  equations  as  is  described  in  reference  5 for  the  plate  vibration  problem.  This 
procedure,  slightly  modified,  is  used  here  for  the  buckling  of  rectangular  cylinders. 

The  energy  principle  from  which  the  governing  equation  (2)  can  be  derived  is  the 
minimization  of  the  following  integral. 


5W.C.  Walton,  Jr.;  Application  of  a General  Finite  Difference  Method  for  Calculating 
Bending  Deformations  of  Solid  Plates;  NASA  TN  D-536;  1960. 


The  origin  of  the  coordinates  yj  and  thus  t?j  has  been  shifted  from  the  center  of  the 
sides  to  the  corner  of  the  cylinder  for  this  finite  difference  analysis.  The  geometrical 
vmr'  atry  of  the  rectangular  cylinder  makes  it  possible  to  treat  only  or  ■ quarter  of  the 
cylinder,  as  shown  in  Figure  3,  where  the  coordinate  system  to  be  used  is  also  shown. 
Here,  as  in  the  previous  analysis,  initially  the  side'  are  treated  as  though  they  are  not 
attached;  thus  the  energy  is  written  in  nondimensional  form  as: 


where  e3  = 0 1 2 /D  n and  e4  = D3  3 /D, , and  the  other  parameters  are  as  previously  defined. 
The  assumption  of  simply  supported  ends  at  x=0  and  x=a  allows  specification  of  w,  and 
w2  as 


w,(£,  T?i)  = w, (77 ! )Sin(n tr£) 
w2(£,  t}2)  = w2  (tj2  )Sin(n7rf) 


(15) 


Subsitution  of  these  expressions  into  the  energy  and  carrying  out  the  integration  with 
respect  to  £ transforms  (14)  to  the  following  one-dimensional  energy  expression 


b/a 


-c/a 


— Ed?7 1 + 

2 


where 
E = 


njr 

T 


4 (»,)■  -2 


, Nlntr)2  2 
+ ' 


(16) 
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2Ay,  e4(n7r)2 
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dt?? 

2 

dr?,  I 

(17) 
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over  0 < 77,  < b/a  with  a similar  expression  over  0 < t?2  < c/a  and  with  w,  and  t?, 
replaced  respectively  by  w2  and  t?2. 

Reduction  to  Algebraic  Eigenvalue  Problem: 

To  carry  out  the  reduction  of  equation  (16)  to  an  algebraic  eigenvalue  problem  using 
finite  difference  procedures,  a grid  is  established  as  shown  in  Figure  3.  In  this  figure 
grid  points  on  the  structure  are  denoted  by  open  symbols,  and  the  grid  points  not  on 
the  structure  required  for  evaluation  of  derivatives  at  the  boundary  points  are  denoted 
by  solid  symbols. 

On  each  side  m grid  points  are  used;  these  are  numbered  2 through  (m+1)  on 

0 < y,  < b and  (m+3)  through  (2m+2)  on  0 < y2  < C.  The  four  off-structure  points 
are  numbered  1,  (m+2),  (2m+3),  (2m+4).  The  first  step  in  the  reduction  is  the  evaluation 
of  the  integrals  in  equation  (16)  in  terms  of  the  energy  density  function  E evaluated 
at  the  grid  points.  This  is  accomplished  using  Simpson's  Rule  for  integration  and  the 
expression  for  the  energy  (16)  becomes: 

1 (|E2  + 4E,  + 26,  + 4ES  + ■ • • 2Em.,  ♦ 4Em  + E^,  ] Al 

2 1 3 

+ (Em+3  + 4Em+4  + 2Em+5  + 4Em+6  + ' ' * 2E  2m  + 4E2m+1  + E2m+2  j ^ (18) 


Ai?j  is  the  grid  point  spacing. 


To  complete  the  reduction  to  an  algebraic  eigenvalue  problem  the  energy  density  function 
Ej  at  any  point  is  expressed  in  terms  of  the  grid  point  displacements  w,-  by  using  finite 
difference  approximations  for  the  derivatives  in  equation  (17)  to  give: 


E:  = (rm)4  (wj)2  - 2e3(n7r)2  (wf) 


wi.1-2wi4wi^1 

Ar?i 


+ e2 


At?? 


+ 4e4  (n  tt)  2 


(Wj.l+Wi+l)2 

4 At?  2 


+ N(nr)2w? 


for  i = 2 to  (m+1)  and  a similar  expression  for  i = (m+3)  to  (2m+2)  with  At?,  replaced 
with  At?2.  The  quadratic  expression  (19)  can  be  written  in  matrix  form  as 
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where  aj  j=1,6  and  dj  are  given  in  Appendix  il.  Substitution  of  these  expressions  into 
equation  (18)  yields  the  following  matrix  expression  for  the  total  energy: 

(w  J 1 [S]  ( w]  + N (w  ) T [T]  [ w j (21) 

w is  a vector  whose  elements  are  the  grid  point  displacements  wj  and  [S]  and  [T]  are 
matrices  constructed  from  the  3 X 3 matrices  in  equation  (20).  The  elements  of  these 
3X3  matrices  are  given  in  Appendix  II. 

Application  of  Constraints: 

Before  the  energy  (21)  can  be  minimized  to  yield  the  equilibrium  equations,  the 
constraints  which  express  the  symmetry  conditions  and  the  conditions  which  join  the  sides 
of  the  cylinder  must  be  imposed.  This  is  accomplished  through  a direct  application  of 
the  procedure  described  in  reference  4. 

From  the  results  of  the  direct  solution  it  was  known  that  the  lowest  buckling 
eigenvalue  for  this  problem  is  associated  with  a deformation  shape  in  which  all  sides  undergo 
symmetric  deformation.  Thus  the  symmetry  conditions  at  i7,=b/a  and  Tj2=c/a  are  the 
vanishing  of  the  slopes  at  these  points.  These  constraints  can  be  expressed  mathematically 
in  terms  of  the  finite  difference  approximations  as 

w2m+3  ~ wm  = 0 

w2m+4  ~ w2m+1  = 0 (22) 

The  constraints  required  to  impose  the  joining  of  the  sides  are  the  vanishing  of  the  normal 
displacements  at  171  =0  and  r72  =0  and  the  continuity  of  the  slope  and  curvature  of  the 
corner.  These  are  expressed  in  finite  difference  form  as: 

w2  =0 

wm+3=0 

1 /2 A77 , (w3-Wi ) - 1/2A772(wm+4— wm+2)  = 0 

1/Atji  (wj  — 2w2+w3)  + 1/Ar?l(wm+2-2wm+3+wm+4)  (23) 
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These  six  homogenous  equations  (22)  and  (23)  can  be  expressed  in  matrix  form  as 

[C]  jw j = 0 (24) 

where  [C]  is  a 6X(2m+4)  matrix  whose  non-vanishing  elements  are  given  in  Appendix  III. 
The  procedure  in  reference  4 provides  a transformation  of  the  vector  |w|  into  a new 
vector  of  generalized  coordinates  [qj  which  satisfy  equation  (24)  and  are  expressed  in 
the  form 


w = [0]  {q  j (25) 

This  transformation  is  used  in  equation  (21)  to  give  an  expression  for  the  energy  in  terms 
of  generalized  coordinates  which  obey  the  symmetry  and  joining  constraints 

{ q )T  IS]  ( q]  + N [q  jT  [T]  jq  j (26) 


where 


[5]  = [0]  T [S]  [0] 

[T]  - m T [T]  [0]  (27) 

Minimization  of  the  energy  with  respect  to  the  generalized  coordinates  jqj  yields 

the  following  equilibrium  equations: 

(S)  (qj  + N (T]  (q  j = 0 (28) 

Solutions  to  these  equations  were  obtained  by  the  use  of  the  threshold  Jacobi  Numerical 
Technique.  Once  the  eigenvalue  N and  the  corresponding  eigenvectors  IQ]  are  obtained, 
the  eigenvectors  can  be  transformed  back  to  the  physical  coordinates  jwj  by  the  use  of 
equation  (25). 

Treatment  of  Nonuniform  Load: 

To  this  point  the  development  of  the  finite  difference  solution  has  considered  the 
uniformly  loaded  cylinder.  The  geometry  of  the  rectangular  cylinder  suggested  that  the 
treatment  of  the  nonuniform  loads  might  be  required. 

A ncnuniform  load  implies 'that  the  parameter  N in  equation  (14)  is  not  a constant 
but  a function  of  tj,  and  r?2.  This  has  the  principal  effect  of  modifying  the  matrix 
element  dj  in  equation  (20)  and  the  interpretation  of  the  buckling  eigenvalue  N.  The 
types  of  load  distribution  of  interest  are  those  in  which  a greater  portion  of  the  load 
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is  supported  at  the  comers  as  illustrated  in  Figure  4.  Such  a distribution  is  conveniently 
modeled  by  two  second  order  polynominals.  In  choosing  the  exact  form  for  these 
polynomials,  it  is  desirable  to  place  certain  constraints  on  their  behavior,  namely,  they 
must  be  continuous  at  the  corner,  take  on  their  minima  at  r?,=b/a  and  r?2=c/a  and  be 
symmetrical  about  these  points.  Functions  satisfying  these  constraints  are: 

N = N ^1  + b,  l^i-b/a)2  - 1/3(b/a)J]j  on  (X?h  <fc/a 

N = N |l  +b,  [{T7,  —c/a)2  + 1/3(c/a)J]|  onO<T?2<c/a  (29) 

The  parameter  FT,  v/hid  vii!  be  the  eigenvalue  found  by  the  solution,  is  the  average 
value  of  the  buckling  st’ess;  l « its  product  with  the  perimeter  gives  the  total  load  applied 
to  the  cylinder  as  is  the  case  v».th  the  eigenvalue  associated  with  the  uniform  load.  This 
facilitates  comparison  of  the  results  from  the  two  loading  conditions,  since  they  are  both 
proportional  to  the  total  load-cariying  capability  of  the  cylinder.  The  factor  b!  appearing 
in  the  expression  for  the  buckling  stress  (29)  is  specified  so  that  the  ratio  of  the  maximum 
value  of  N to  its  minimum  value  takes  on  the  desired  value. 

The  changes  in  the  functional  form  of  the  buckling  stress  distribution  described  above 
changes  the  element  of  the  matrix  [T]  in  equation  (21),  or  more  specifically,  it  changes 
the  formula  used  to  compute  dj  in  equation  (20).  The  formulas  to  be  used  in  this  case 
are  presented  in  Appendix  II.  The  form  of  the  equation  to  be  solved  (28)  remains 
unchanged;  however,  the  elements  of  (T]  are  changed  and  the  interpretation  of  the 
eigenvalue  R becomes  that  of  the  average  buckling  stress. 
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Figure  4 Illustration  of  the  nonuniform 
load  distribution 
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EXPERIMENTAL  DETERMINATION  OF  THE  BUCKLING  LOAD 


The  validity  and  range  of  application  of  the  theoretical  solutions  presented  in  the 
previous  sections  can  be  established  by  correlation  with  experimental  results.  To  make 
this  possible,  a series  of  rectangular  fiberboard  cylinders  were  tested  to  produce  a body 
of  experimental  data  with  which  the  theoretical  predictions  could  be  compared.  The 
experimental  procedures  used  are  descr1-  here. 

Each  specimen  was  fabricated  in  the  form  of  a rectangular  cylinder;  that  is,  a 
rectangular  container  with  no  top  or  bottom  flaps,  from  V3C  grade  corrugated  fiberboard 
with  the  corrugations  running  parallel  to  the  cylinder  axis.  Each  of  the  cylinders  had 
a stapled  lap  seam  in  one  corner  parallel  to  the  axis  of  the  cylinder.  The  sizes  of  the 
cylinders  can  be  expressed  by  the  two  aspect  ratios  b/a  and  c/a  and  magnitude  of  the 
cylinder  height  a.  The  aspect  ratios  of  the  specimens  tested  ranged  from  0.2  to  0.8  with 
cylinder  heights  of  0.457,  0.609  and  0.762  m included. 

The  buckling  tests  were  carried  out  on  a Tinius  Olsen  testing  machine  capable  of 
accepting  the  large  specimens  involved  in  this  testing  program.  A photograph  of  this 
testing  machine  with  a specimen  in  position  is  shown  in  Figure  5. 

The  tests  were  conducted  by  loading  the  specimen  starting  at  the  unloaded  state 
and  increasing  the  load  at  a uniform  rate  of  1.0  cm/min.  Loading  was  continued  until 
deformation  of  the  specimen  occurred  w:th  no  further  increase  in  load.  The  buckling 
load  was  taken  to  be  the  maximum  load  reached  as  measured  on  the  Tinius  Olsen 
load-measuring  system.  The  total  load  was  reduced  to  a buckling  stress  by  dividing  this 
buckling  load  by  the  perimeter  of  the  cylinder.  For  each  size  cylinder  two  test?  were 
performed. 
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RESULTS 


The  presentation  of  the  results  will  begin  with  a discussion  of  the  experimental 
behavior  of  the  buckling  stress  as  a function  of  cylinder  geometry  and  follow  with  a 
discussion  of  the  theoretical  predictions  and  a comparison  of  these  predictions  with  the 
experimental  behavior.  The  section  will  include  the  discussion  of  the  validity  of  the  theory 
and  some  suggestions  for  its  modification. 

Experimental  Results 

The  graphs  in  Figures  6 and  7 illustrate  the  behavior  of  h e buckling  stress  as  a 
function  of  aspect  ratio  and  cylinder  height.  In  both  figures  the  nondimensional  buckling 
stress  is  plotted  as  a function  of  aspect  ratio,  b/a.  The  buckling  stress  is  nondimensionalized 
according  to  equation  (2)  with  D,  j = 14.77  N-m,  as  given  in  reference  6.  In  Figure  6 
results  are  given  for  two  values  of  the  other  aspect  ratio,  c/a,  with  the  cylinder  height 
constant;  while  in  Figure  7 results  are  given  for  three  values  of  cylinder  height  with  the 
aspect  ratio,  c/a,  held  constant.  The  results  in  Figure  6 show  that  the  buckling  stiess 
increases  with  decreases  in  either  of  the  aspect  ratios.  Thus,  for  buckling  failure  a 
rectangular  cylinder  becomes  stronger  as  the  width  of  the  sides  become  smaller  in 
comparison  to  the  height. 

A similar  behavior  with  respect  to  uipect  ratio  b/a  is  seen  from  the  results  in  Figure  7. 
Additionally,  it  is  seen  that  the  buckling  stress  decreases  with  the  cylinder  height.  This 
is  an  unexpected  result.  The  nondimensional  governing  equations  and  boundary  conditions 
contain  the  nondimensional  buckling  stress  and  the  aspect  ratios,  thus  since  the  cylinder 
height,  a,  does  not  appear  explicitly  one  would  expect  the  nondimensional  buckling  stress 
to  be  independent  of  cylinder  height.  That  is,  that  the  nondimensional  buckling  stress 
be  completely  specified  by  the  aspect  ratios  c/a  and  b/a  for  all  values  of  a.  As  a result 
one  expects  the  total  dependence  of  the  physical  buckling  stress,  Nx,  on  the  cylinder 
sizes  as  measured  by  thr.  cylinder  height,  a,  to  be  expressed  in  the  definition  of  the 
nondimensional  buckling  stress: 


0.. 


Since  one  would  expect  that  N would  be  independent  of  "a"  the  experimental  behavior 
illustrated  in  Figure  7 then  suggests  some  inadequacy  in  the  theory  which  will  be  discussed 
subsequently  in  connection  with  the  comparison  of  theoretical  and  experimental  results. 


6E.C.  Steeves;  Orthotropic  Bending  Properties  of  Fiberboard,  Technical  Report 
75-97-AMEL;  U.S.  Army  Natick  Development  Center;  1975. 
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Figure  7 Experimental  buckling  stress  as  a function  of 
aspect  ratio  and  cylinder  height 
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Theoretical  Results  and  Comparison  with  Experiment 

It  will  be  recalled  that  solutions  to  the  rectangular  cylinder  buckling  problems  were 
obtained  for  a uniform  load  by  both  exact  and  approximate  techniques  and  for  a 
nonuniform  load  by  approximate  techniques.  The  results  of  these  solutions  are  shown 
in  Figures  8 and  9,  as  plots  of  the  nondimensional  buckling  load  as  a function  of  the 
aspect  ratio.  These  results  were  obtained  using  the  bending  stiffness  for  V3C  fiberboard 
given  in  reference  6.  In  Figure  8,  which  gives  results  for  cylinders  having  the  aspect 
ratio  c/a  * 0.6  with  Va  as  a variable  ranging  from  0.2  to  0.8,  a comparison  is  shown 
of  the  exact  solutions  and  the  approximate  solutions  for  the  buckling  stress.  These 
solifions  are  everywhere  vrthin  ri % of  each  other,  which  is  sufficiently  lose  to  establish 
the  validity  of  the  approximate  solution. 

The  experimental  results  are  also  shown  on  Fiqure  8.  Comparison  of  these  results 
with  the  theoretical  predictions  based  on  a uniform  load  indicates  poor  agreement.  The 
theory  has  the  proper  trend  of  decreasing  buckling  stress  with  increasing  aspect  ratio, 
but  the  theory  underestimates  the  buckling  stress  by  as  much  as  a factor  of  two.  Thus 
the  theoretical  model  using  a uniform  load  does  not  give  satisfactory  results. 

The  results  presented  in  reference  2 indicate  the  existence  of  a nonuniform  load 
during  buckling,  and  this  seems  reasonable  when  one  considers  the  relatively  greater  stiffness 
of  corner  regions  of  the  rectangular  cylinder  in  comparison  to  that  of  the  central  regions 
of  the  sides.  Being  stiffer,  the  corners  will  suoport  a greater  portion  of  the  load  when 
the  cylinder  is  subjected  to  a uniform  compression  or  axial  shortening  which  is  the  type 
of  loading  that  occurred  during  the  buikling  tests.  Because  of  this,  solutions  corresponding 
to  the  experimental  results  were  obtained  using  the  nonuniform  load  distribution  specified 
in  equation  (29).  It  was  not  possible  to  obtain  Jhese  solutions  in  a direct  manner  because 
this  load  distribution  contains  two  unknowns,  N,  the  average  value  of  the  buckling  stress, 
and  b|,  j parameter  which  determines  the  ratio  of  the  maximum  stress  to  the  minimum 
stress.  N will  be  determined  as  the  eigenvalue  of  the  solution  but  b,  must  be  specified. 
For  the  results  presented  in  Figures  8 and  9,  b|  was  chosen  so  that  the  results  were 
in  fair  agreement  for  one  value  of  b/a  and  the  same  value  of  b)  used  for  the  remaining 
values  of  b/a.  The  values  of  b,  used  are  giver,  on  the  figures  and  are  different  for  the 
two  values  of  the  aspect  'atio  c/a  Examination  of  the  results  presented  in  Figure  8 
show  that  the  predictions  of  the  average  buckling  stress  based  on  the  described  load 
distribution  are  in  fair  agreement  with  experimentally  determined  values  of  the  buckling 
stress.  The  experimental  values  can  also  be  considered  as  average  values  because  they 
were  computed  by  dividing  the  total  buckling  load  by  the  perimeter  of  the  cylinder.  These 
comparisons  between  theory  and  experiment  indicate  that  behavior  can  be  characterized 
with  a nonuniform  load  distribution.  However,  the  exact  form  of  this  distribution  is 
not  known  and  cannot  in  an  exact  way  be  determined  from  the  procedure  developed 
here.  This  suggests  treatment  of  the  geometrically  nonlinear  buckling  problem  in  which 
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the  rectangular  cylinder  is  subject  to  a uniform  compression.  The  buckling  stress 
distribution  would  be  one  of  the  results  of  this  problem  instead  of  a specified  input. 

This  nonuniform  load  result  also  suggests  an  explanation  of  the  unexpected  variation 
of  the  nondimensional  buckling  stress  with  cylinder  height  found  experimentally  and 
discussed  previously.  The  distribution  of  load,  and  thus  stress,  resulting  from  a uniform 
compressive  deformation  is  dependent  on  the  distribution  of  stiffness  along  perimeter  of 
the  cylinder.  As  indicated  above,  this  distribution  is  influenced  by  the  presence  of  the 
corners  but  it  is  also  believed  to  be  influenced  by  the  height  of  the  cylinder.  For  example; 
if  compressive  behavior  predominates  in  the  region  near  the  corner  and  bending  behavior 
in  the  central  region  of  the  sides,  then  the  corner  stiffness  will  vary  inversely  with  cylinder 
height  and  the  central  stiffness  inversely  with  the  square  of  the  cylinder  height.  As  a 
result,  as  the  height  increases  the  stiffness  of  the  comer  will  become  more  predominant 
resulting  in  a concentration  of  load  at  the  corner  with  the  resulting  increase  in  the  buckling 
strength  with  cylinder  height.  This  mechanism  provides  behavior  in  agreement  with 
experimental  results  and  illustrates  a possible  explanation  for  the  variation  of  buckling 
stress  with  cylinder  height.  It  is  beIJeved  that  the  true  mechanism  behind  this  behavior 
will  be  revealed  by  solution  of  the  nonlinear  buckling  problem  which  includes  the 
interaction  of  bending  and  compressive  modes  of  deformat:on. 

Figure  9 illustrates  the  behavior  of  the  buckling  stress  as  a function  of  the  aspect 
ratio,  b/a,  with  the  other  aspect  ratio,  c/a,  equal  to  0.3.  The  theoretical  predictions  from 
both  the  uniform  and  nonuniform  load  are  considerably  different  in  relation  to  the 
experimental  results  from  that  found  for  c/a  = 0.6.  For  the  present  case,  c/a  = 0.3, 
the  theoretical  nonuniform  load  results  are  in  fair  agreement  with  the  experimental  results 
for  0.6  < b/a  < 0.8.  For  b/a  < 0.6  the  theoretical  results  are  in  poor  agreement  with 
large  overprediction  of  the  buckling  stress.  The  phenomenon  behind  this  behavior  is 
suggested  by  the  buckling  behavior  of  flat  plates  subjected  to  stress  beyond  the  elastic 
limit  shown  as  an  insert  in  Figure  9 and  taken  from  reference  7.  This  data  shows  that 
as  long  as  the  buckling  stress  is  below  the  yield  stress,  the  classical  buckling  result  holds; 
once  the  yield  stress  is  reached,  the  plate  will  buckle  regardless  of  size.  The  results  for 
the  rectangular  cylinder  shown  on  Figure  9 exhibit  a similar  behavior,  taking  the  solution 
under  nonuniform  load  as  the  classical  buckling  result  and  the  experimental  results  as 
the  branch  representing  the  stability  boundary  due  to  plastic  buckling.  In  this  case  buckling 
at  a constant  stress  level  is  not  observed  because  fiberboard  does  not  have  a sharply  defined 
yield  point,  and  the  state  of  stress  is  not  uniform  as  is  the  case  in  the  example  from 
reference  7.  To  further  confirm  the  existence  of  plastic  buckling  in  the  experimental 
results,  a measure  of  the  level  of  the  average  stress,  N,  at  which  the  stress  in  the  rectangular 
cylinder  reaches  the  yield  point  is  needed.  A convenient  measure  is  provided  by  the 
maximum  stress  occurring  in  the  nonuniform  distribution  used  in  this  analysis.  The  value 
of  N for  which  this  maximum  stress  reached  the  yield  point  was  computed  and  is  plotted 
on  Figure  9.  As  can  be  seen  this  value  is  in  good  agreement  with  the  value  at  which 

7S.P.  Timoshenko  and  JM.  Gere;  Therory  of  Elastic  Stability;  McGcaw-Hil1  Book  Gompany, 

1961. 
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the  experimental  and  theoretical  results  depart.  Since  this  value  of  N is  only  a measure 
of  fhe  point  at  which  plastic  behavior  begins,  it  is  not  surprising  that  the  expe'imental 
results  are  greater.  Thus  it  appears  that  in  addition  to  the  need  for  treatment  of  the 
geometrically  nonlinear  buckling  problem,  there  is  also  a need  to  consider  inelastic  buckling 
in  the  case  of  rectangular  cylinders  with  small  aspect  ratios. 

These  results  have  substantial  impact  on  the  selection  and  specification  of  fiberboard 
grades  for  a given  packaging  situation.  This  impact  comes  from  the  fact  that  for  the 
buckling  mode  of  failure  the  bending  stiffness  is  the  materiel  property  that  determines 
package  performance.  That  is,  an  increase  in  the  bending  stiffness  gives  a proportional 
.r, cruise  in  the  buckling  sti ongth  as  shown  by  equation  (2d).  The  p^^ent  design  and 
specification  parameter  for  fiberboard  is  the  Mullen  burst  strength,  a property  guide 
different  from  the  bending  stiffness.  Thus  it  is  recommended  that  the  bending  stiffness 
of  fiberboard  be  used  as  an  additional  specification  parameter.  Such  a parameter  could 
be  determined  by  simplification  of  the  bending  test  procedure  described  in  reference  6. 
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Figure  8 Theoretical  predictions  of  the  buckling 
stress  and  comparison  with  experiment 


29 


b/a, ASPECT  RATIO 

Figure  9 Theoretical  prediction  of  the  buckling 
stress  and  correlation  with  experiment 
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CONCLUSIONS 


The  results  of  a theoretical  and  experimental  investigation  of  the  buckling  of 
rectangular  cylinders  has  been  presented.  As  a result  of  the  comparison  of  the  theoretical 
results  with  those  determined  experimentally,  the  following  conclusions  can  be  reached: 

1.  A theoretical  model  based  on  a uniform  loading  situation  is  inadequate  for 
predicting  the  buckling  stress  of  rectangular  cylinders. 

2.  The  classical  buckling  theory  using  a nonuniform  load  shows  promising 
correlation  with  experimental  results;  however,  with  the  classical  theory  it  is  not  possible 
to  completely  determine  the  load  distribution.  To  remedy  this,  the  geometrically  nonlinear 
buckling  with  a uniform  compression  or  shortening  of  the  cylinder  should  be  studied. 

3.  For  rectangular  cylinders  w.th  small  aspect  ratios,  inelastic  buckling  or  material 
nonlinearity  must  be  Included  in  the  theory  to  give  adequate  prediction  of  the  buckling 
strength. 

4.  The  bending  stiffness  should  be  used  as  an  additional  design  and  specification 
parameter  for  fiberboard  used  in  packages  and  containers. 
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APPENDIX  I 


None  vanishing  elements  of  the  matrix  U 


U , [ = 0Cos0b/a 
Ul2  = 0Sin0b/a 
U,  3 = -|3Cos0c/a 
Ul4  = 0Sin0c/a 
U2l  = -02Sin|3b/a 
U22  = /32Cos0b/a 
U2  3 = -/32  Sin/3c/a 
U24  = -02Cos0c/a 
U33  = 0Cos0c/a 
U 3 4 = 0Sin0c/a 
U35  = -/3Cos0b/a 
U36  = 0Sin/3b/a 
U4  3 = */32  Sin/3c/a 
U44  = 02Cos0c/a 
U45  = -02Sin0b/a 
U46  = -02Cos0b/a 
U55  = 0Cos/3b/a 
U56  = 0Sin/3b/a 


+ f 2 aCoshab/a 

- fiOSinhob/a 

- f4aCoshac/a 

- f3oSinhcrc/a 
+ f2a2Sinhab/a 

- f,o2Coshab/a 
+ f4a2Sinhac/a 
+ f3o2Coshocc/a 
+ f4aCoshac/a 

- f3c*3inhac/a 

- f2  aCoshab/a 

- f|OSinhab/a 
+ f4o2Sinhac/a 

- f3a2Coshac/a 

+ f2a2Sinhab/a 
+ fiOJCoshab/a 
+ f2  aCoshab/a 

- f,  ctSinhab/a 
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APPENDIX  I 
(cont'd) 


U5  7 = -0Cos0c/a 
U5g  = /JSin(3c/a 
U65  = -(32Sin/3b/a 
U66  = /32  Cos^b/a 
U67  = -02Sin0c/a 
U68  = -02Cos0c/a 
U7 , = -0Cos/3b/a 
U72  = 0Sinj3b/a 
U17  = 0Cos|3c/a 
U78  = 0Sin0c/a 
U8 1 = -02Sin/3b/a 
U82  = (32  Cos0b/a 
U87  = -02Sinj3c/a 
U8  8 = (32Cos/3c/a 


- f4orCoshac/a 

- f3aSinhac/a 
+ f2a2Sinhab/a 

- f,a2  Coshab/a 
+ f4a2Sinhac/a 
+ f 3a2  Coshac/a 

- f2aCoshab/a 

- f,aSinhab/a 
+ f4aCoshac/a 

- f3aSinhac/a 
+ f2aJSinhab/a 
f f,a2Coshab/a 
+ f4a:Sinhac/a 
+ f3  a2  Coshac/a 


APPENDIX  II 


a i 
a2 
a3 
a4 

35 

a6 

di 

At?, 


1/2(e2/AT?*  + e4(njr)2/At??J 
1/2  ^3(njr)2/AT?^  - 2e2/Ar,jj 
1/2[e2/Ar?t  - e4(nTr}2/Ai?j] 

1/2  |nff)4  + 4e3  (nw)2  /Ai?^  + 4e2/At?f] 
1/2  (n?r>2 /At?2  - 2e2/Ar?t] 

1 /2  ^e2  /At?“?  + e4(nTr>2/AT?,J 
1/2  (nff)2  uniform  buckling  stress 
At?2  where  applicable 


1 + 


(c/a)2  b 1 b b 

b'  _ ->2J  " “I  )21  0<n,  <- 

(b/a)2  a 3 a a 


di  = < 


1 + 


b,  ((T?2j--)2  - - A2) 

a 3 a 


0<t?2  <1 
a 


nonuniform  buckling  stress 


? 
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APPENDIX  III 


Non-vanishing  element  of  the-  constraint  matrix. 
C(1,2)  = 1.0 
C(2,  m+3)  = 1.0 
C(3,3)  = 1/2  At,, 

C(3,1)  = -1/2AIJ, 

C(3.  m+4)  = -1/2Atj2 
C(3,  m+2)  = 1/2Atj, 

C(4,1)  = 1 /At?  * 

C(4,2)  = -2/At,] 

C(4,3)  = 1/At,? 

C(4,  m+2)  = 1/At,? 

C(4,  m+3)  = -2/At,? 

C(4,  m+4)  = 1/At,? 

C(5,  2m+3)  = 1 
C(5,  m)  = -1 
C(6,  2m+4)  = 1 
C(6,  2m+1)  = -1 
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LIST  OF  SYMBOLS 


a 


Ai.Bj.Cj.Dj 

b,c 

b, 

dj 

Djk 


ei  ,e2 ,e3,e4 
E 

*i 

FjlTJj) 

m 

n 

Nx 

N 

R 

p 


Wj 


x.Vj 


a,(3 


Cylinder  height 
Stiffness  matrix  constants 
Constants  of  integration 
Cylinder  side  widths 

A constant  used  to  describe  the  nonuniform  buckling  load 

Buckling  load  matrix  constant 

Bending  stiffness  constants 

Bending  stiffness  ratios 

Energy  density  function 

Consents  defined  by  equation  10 

Function  describing  the  normal  displacement 

Number  of  grid  points  per  side 

Axial  half-wave  number 

Axial  stress  resultant 

Nondimensional  axial  stress  resultant 

Average  value  of  N 

Product  of  eigenvalues 

Normal  displacement  of  the  ith  side  of  cylinder 
Position  coordinates  on  the  ith  side  of  cylinder 
Characteristic  numbers  defined  by  equation  6 
Nondimensional  coordinates  corresponding  with  (x,yj) 
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LIST  OF  SYMBOLS  (corn'd) 


Matrices 

Constraint  matrix 
Matrix  of  eigenvectors  jq] 

Stiffness  matrix  associated  with  |wj 
Stiffness  mat  ix  associated  with  jq  J 
Buckling  loaa  matrix  associated  with  jwj 
Buckling  load  matrix  associated  with  [qj 
Boundary  condition  matrix 
Matrix  defining  the  generalized  coordinates 
Vectors 

Vector  of  integration  constants 

Vector  of  finite  difference  grid  point  displacements 

Vector  of  generalized  coordinates  defined  by  equation  , 5) 
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